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A relativistic light front formulation of nuclear dynamics is applied to infinite nuclear matter. A 
hadronic meson-baryon Lagrangian, consistent with chiral symmetry, leads to a nuclear eigenvalue 
problem which is solved, including nucleon-nucleon NN correlations, in the one-boson-exchange 
approximation for the NN potential. The nuclear matter saturation properties are reasonably well 
reproduced, with a compression modulus of 180 MeV. We find that there are sa 0.05 excess pions 
per nucleon. 
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An important class of high-energy experiments involv- 
ing nuclear targets is most easily understood using light- 
front variables. Consider the lepton-nucleus deep inelas- 
tic scattering experiments [H which found a significant 
difference between the parton distributions of free nu- 
cleons and nucleons in a nucleus. This is interpreted 
as a ~ 10% shift in the momentum distribution of va- 
lence quarks towards smaller values of the Bjorken vari- 
able xbj- The Bjorken variable is a ratio of the plus- 
momentum p + = p° + p 3 of a quark to that of the tar- 
get. Treating the nucleus as a collection of nucleons gives 
x Bj = P + /k + , where k + is the plus momentum of a 
bound nucleon. The spatial variable canonical to k + is 
x~~ = x° — x 3 , and the time variable is x° + x 3 = x + . To 
use these variables is to use the light front variables of 
Dirac 0. See the recent reviews ||. 

Deep inelastic scattering depends on the probability 
f{k + ) that a bound nucleon has a momentum k + . Other 
nuclear reactions, such as (e,e') and (p,2p) depend also 
on this very same probability Q. The quantity f{k + ) 
is simply related to the square of the ground state wave 
function, provided light front dynamics are used. In the 
equal-time formulation f(k + ) is obtained as a response 
function that depends on matrix elements between the 
ground and all excited states. 

The use of light front variables is convenient for in- 
terpreting experiments, but one is faced with computing 
the ground state nuclear wave function using light front 
variables. We introduce here a theory of infinite nuclear 
matter, using light front dynamics, in which the effects of 
correlations are taken into account. This is a light front 
Brueckner theory of nuclear matter. 

Light front dynamics is a three-dimensional Hamilto- 
nian formalism. A Lagrangian is chosen and one derives 
field equations which allow the elimination of dependent 
degrees of freedom. The light front Hamiltonian P~ is 
obtained from (J 1 ^ is the energy momentum ten- 

sor), and P + , the plus- momentum operator is obtained 
from T ++ . We use a chiral model, detailed in Refs. H|), 
in which the nuclear constituents are nucleons ip, pions, 
scalar mesons, vector mesons, and 77 mesons. 

We start by briefly reviewing some of the necessary for- 
mal details presented in Ref. B (which presented a mean 
field calculation of the properties of nuclear matter) . Our 
specific equations employ the rest frame of nuclear mat- 
ter. The derivation of is accomplished by solving 
for the dependent nucleon fields in terms of the indepen- 
dent ones. The simple completion of this task requires 
setting the + components of the vector fields V 11 to zero. 
This is accomplished using a transformation |7) on the 
nucleon fields, in which one replaces V^ L with other fields 
(y+ — 0) which enter in the nucleon field equations. 
We carry this feature through in our calculations (includ- 
ing the related mean fields), but drop the bar to simplify 
the notation. More detailed explanations are given in 
Refs. [§ and @. 



We now confront the eigenvalue problem for the ground 
state of infinite nuclear matter: 

P-\V) = M A \y). (1) 

For a nuclear system at rest we must also satisfy P + Y$>) = 
Ma\^). The operator P~ can be divided into kinematic 
P ~(N) and interaction J terms: 

P-=P (N) + J (2) 

in which P ~(N) is the kinetic contribution to the P~ 

operator, giving p± p™ for the minus-momentum of free 
nucleons. The operator J is the sum of three terms J = 
i>i + i>2 + ^3 : v i gives the single meson-nucleon vertex 
functions; 1*2 accounts for meson emission followed by 
instantaneous nucleon propagation followed by another 
meson emission; and, v% accounts for the instantaneous 
propagation of vector mesons. 

The traditional path of using a two-nucleon potential 
and temporarily eliminating the meson degrees of free- 
dom by adding and subtracting a two-nucleon potential 
to the Hamiltonian is followed. The terms involving the 
difference between J and the two-nucleon potential are 
handled later using perturbation theory. Then the ef- 
fective Lagrangian is £/c = V>(?7 • d — M)ip — ^,with the 
two-nucleon interaction density /C/2. Using gives the 
corresponding P ± operators, with PZ = P Q ~(N) + K. K 
is the volume integral of /C, and 

P- =P K + J-K + Po(m), (3) 

where P ~(m) accounts for the non-interacting mesonic 
contribution to P~ . 

The purely nucleonic part of the full wave function is 
defined as | $), and is the solution of the light front 
Schroedinger equation 

Pk I $ ) = ( P o W +K)\$) = M I $). (4) 

The interaction K is obtained using second-order per- 
turbation theory. Consider the scattering process 1 + 2 — > 
3 + 4. We find 

(3,4|Jf|l,2> - (3,41^(^7)^ +«3|1, 2), (5) 
with go(P^j) = P -\p- i where P~j is the negative com- 

ij 

ponent of the total initial (12) or final (34) momentum 
which here are the same. 

We discuss the energy denominator (P^ — P^)k + . 
Suppose that > k^ . Then the emitted meson of mass 
/i and momentum k, with k + = k^ — k^ , k± = ki^—k3j_, 

k- = and k+(Pr o - P ~) = fc+(fcf - fc 3 ") - 

kj 2 + /i 2 = q 2 — /1 2 , where q = k\ — k%. That this form is 
the Lorentz-invariant inverse of the Klein-Gordon propa- 
gator indicates that K accounts for the nucleon-nucleon 
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one-boson exchange potential. Note that q 
with q 2 accounting for retardation effects. 

The interaction K is strong, so that we sum the uni- 
tary set of iterations to obtain the two-nucleon integral 
equation for the transition matrix T: 

(3,4|T|1,2) = (3,4|V|1,2) + / ^(3,4|7|5,6) 

d 2 p±da 



A5 , Afi 

2M 2 -> 2 » 



a(l — a) pi _ 



p 2 ± +M 2 
ij a(l-a) 



up the determinant \<j> > are to be given in Eq. (9), 
which is not the Dirac equation. The self-consistent fields 
U which enter are to be obtained from the expectation 
value of the G-matrix in | ifi) . The use of symmetries E3] 
causes the only non-vanishing mean-fields to be scalar 
Us and vector fields Uy. Neglect of the very small JTJ] 
spatial component of Uy leaves us with Uy. 

We now obtain the Slater determinant | <j>) and the 
related scalar Us and vector Uy potentials. The eigen- 
(5,6|T|1,2). (6) values and eigenvectors of the resulting light-front mean 
field Hamiltonian P^ F are easy to obtain, and can be 



The plus-momentum variable is expressed in terms of a 
momentum fraction a such that p$ — ceP^, and p± = 
(1 — a)j>5i_ — ap6±- The potential V is a kinematic factor 
times K . 

We mention again that V is derived using the one 
boson exchange approximation, so that effects of multi- 
mesons in flight are not included in Eq. (6). The effects 
of the stretched box diagram has been studied in the 
light front calculations of Ref. || , and can be important 
in maintaining the full Lorentz invariance of the scatter- 
ing amplitude. The present results (6) for the on-shell 
T-matrix are invariant under boosts in the z-direction. 

The Weinberg- type Q equation (||), can be re- 
expressed using the variable transformation ||: a — 
E ^ p f , with E(p) = ^p-p + M 2 . The result is: 

(3,4|T|1,2) = (3,4|7|1,2}+ f ^<3,4|V|5,6) 



A 5 ,A 6 



M 2 d 3 p 
E{p) pf — p 2 + ie 



: (5,6|T|1,2), (7) 



which shows a manifest rotational invariance. 

The one-boson exchange potential V obtained with our 
formalism is constructed elsewhere [§. We choose the 
meson-nucleon coupling constants and include the effects 
of form factors to reproduce the NN data. This light front 
work yields potentials and transition matrix elements es- 
sentially identical to the Bonn OBEP potentials [ p|Jl2| 
because the differences arising from including retardation 
effects, and using Eq. ([?]) are very small. Thus we find 
a good description of the deuteron and the two-nucleon 
scattering data below the inelastic threshold j|] . 

Now that V is determined, we turn to solving Eq. (Q). 
It is worthwhile to begin by explaining how our light 
front quantized calculation of the energy and density of 
nuclear matter includes the effects of two-nucleon NN 
correlations. We employ a light front version of the stan- 
dard technique of obtaining the G-matrix (see Eqs. (11) 
and (12) below) which is a nuclear matter generaliza- 
tion of the T-matrix of Eq. (6) . We then determine the 
nuclear density by finding the Slater determinant \<p > 
which minimizes the expectation value of P~ subject to 
the constraint that the expectation values of P + and P~ 
are equal. The single particle wave functions that make 



MF 

chosen to best approximate the effects of the two-nucleon 
interaction. The mean field light front Hamiltonian den- 
sity T^fp is obtained and its volume integral is a single- 
nucleon operator P MF - ^ one expands ip in terms of 
the single particle wave equation including U as the only 
interaction, one finds that P MF = Pq(N)- The ground 
state eigenvector of this operator is a Slater determinant 
denoted as | 4>): Pq(N) | <j>) — to | tj>), and 



I $> -I 0) + 

where A = 1— | <fi 
are obtained from 
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M -Pu(N)~K 



-AK 



(8) 



. The single nucleon orbitals | i) 



(Pmf - Uy) I »'>+ = (ai ' kix + P(M + U s )) \ 

k+ I »)- = (a ± • fe i± + 0(M + U s )) I i)+, 



P 



MF 



i)- 



(M + Us? 



hi 



i). 



(9) 



The subscripts ± refer to the independent (+) and de- 
pendent (— ) projections of the nucleon fields. 

Using standard manipulations on Eq. (||) , and keeping 
correlations between pairs of nucleons only leads to: 
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(10) 



where Tij is a two-nucleon operator: 
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■Ar tJ K). (ii) 



JV ljl 3 3 Pr-AP (N)A 

The notation i, j refers to the nucleon orbitals of Eq. (||). 
The approximation that U is independent of orbital i 
causes the vector potential to cancel out in the denomi- 



nators e.g. — (pg +Pq). Starting with Eq. ([11]), we 
obtain an equivalent three-dimensional, manifestly rota- 
tionally invariant integral equation by using the medium- 
modified version of the variable transformation |j] : a = 



■2E~ 



,with E* 



pp 



M* 



The result is: 



(3,4|G|1,2) = <3,4|7|1,2)+ f ^(3,4|U|5,6) 
M* 2 d 3 P ' 



A5 ,Ar 



E* p 2 - p 2 



(5,6|G|1,2), (12) 
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where G is related to Y by a kinematic factor. The op- 
erator Q is the two-body version of A and projects the 
momenta p§ and pe above the Fermi sea. The evaluation 
of V using the basis of Eq. is the essential difference 
between relativistic and non-relativistic Brueckner the- 
ory. 

We now specify the Slater determinant | <j>). The oc- 
cupied states are to fill up a Fermi sea, defined in terms 
of a Fermi momentum, kp, that is the magnitude of a 
three vector. This three vector can be defined jl4|] as: 
k + = E*(k) + k 3 , which specifies k 3 . Computing the en- 
ergy and plus momentum then proceeds by taking the ex- 
pectation values of and T ++ . Note that the result- 
ing formulae can be cast in a more familiar form by using 
a discrete representation of the single nucleon states, a 
set of spinors: | a), with a = fc, A such that (a \ a) = 1. 
These light front spinors of Eqs. (9) can be related to 
equal time spinors. Multiply the first of (9) by 7 + and 
the second by j~ . Then add the two equations to obtain 
a new equation. Its solution turns out to be a spinor 
(x | a) which differs from the equal time (ET) spinor by 
only a phase factor, exp {—i Jv ^ ), which cancels out of 
matrix elements under our approximation that Uy is a 
constant. Then the Brueckner light front Hartree-Fock 
(BHF) approximation is to the mean-field is given by 



U{a) = Re(af3 | G \ a(3) a . 

f3<F 



(13) 



The nuclear matter G-matrix is the solution of the in- 



tegral equation, (12) in a plane wave basis. We use the 
implicit definition of k 3 to obtain the total P and rel- 
ative three-momenta of two nucleons interacting in nu- 
clear matter. An angle averaging procedure |l5[ ], is used 
for P and to obtain the Pauli operator Q. The calcu- 
lation of the nuclear matter G matrix involves a self- 
consistency, since the solution of Eq. (|l2] ) for G requires 
knowledge of M * which, in turn, is determined from G 
via Eq. (|l3|). Solutions are obtained using an iterative 
procedure. We obtain the nuclear mass from the average 
1 of P^. M = £ Q<F 6 Q + | E a ,0< F (aP I G | a(3) a . 

The results for the binding energy per nucleon E/A 
in nuclear matter as a function of density, are shown in 
Fig. 1 as the solid line. The curve has a minimum at 
E/A = -14.71 MeV and k F = 1.37 fin -1 , and the in- 
compressibility (compression modulus) is K = 180 MeV 
at the minimum. These values agree well with the em- 
pirical values E/A = -16 ± 1 MeV, k F = 1.35 ± 0.05 
fin -1 , and K = 210 ±30 MeV @. Note also that the ef- 
fect of using non-relativistic Brueckner theory is that the 
saturation density is predicted too high, but the incom- 
pressibility is about the same as ours Jl7j. Comparison 
with the results of an earlier relativistic calculation |Ts| ] 
shows that the only significant difference occurs for the 
incompressibility which is reduced from the earlier value 
of 250 MeV. This is due to using medium-modified single- 



particle energies in the evaluation of the retardation term 
]l9|] , and also to using new values of the mean fields: 
M (kp) = Mpf + Us — 718 MeV and U v = 164.4 MeV. 
These potentials have magnitudes that are considerably 
smaller than those obtained using mean field theory pOl . 
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FIG. 1. Binding Energy per nucleon in nuclear matter, 
E/A as a function of Fermi momentum kp . The solid line 
is our result using light-front Brueckner theory. The dotted 
curve is obtained when the medium effect on meson retarda- 
tion is omitted. The dashed line is the result from nonrel- 
ativistic Brueckner theory. The box represents the range of 
values allowed by extrapolation from data. 

Turn now to the mesonic component of the nuclear 
wave function. The purely nucleonic part of the Fock 
space is | <&) of Eq. (||) , and the full wave function is 
| *> of Eq. qi|). The mesonic content is contained in the 
difference between the two wave functions. Using Eq. (||) 
in Eqs. ([!]) and (Q) allows us to obtain 



*> -I $) 



1 



M A - Aa>P-A, 



-A$(J-X) | $), (14) 



where A$ = 1— | $)($ |. The second term accounts for 
the mesonic content. Taking the expectation value of the 
mesonic number operators yields the average number of 
mesons N m : 



(15) 



where V m represents the potential caused by the ex- 
change of a pion or scalar meson, and its dependence 



on P f . is given in Eq. 



For the pion the result is 



about 5.0% per nucleon [2JJ. This is smaller than the 
18% of Friman et al |^3| because we use scalar mesons 
instead of intermediate A states to provide the bulk of 
the attractive force. 

The values of Us, Uy and allow us to assess the 
deep inelastic scattering of leptons from our version of the 
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ground state of nuclear matter. Using M*(kp) = 744 
MeV Q and neglecting the influence of two-particle- 
two-hole states to approximate f(k + ) shows that nu- 
cleons carry 81% (as opposed to the 65% of mean field 
theory of the nuclear plus momentum. This repre- 
sents a vast improvement in the description of nuclear 
deep inelastic scattering. The minimum value of the ra- 
tio F2A/F2N, obtained from the convolution formula 



F 2A (x) 
A 



dyf(y)F 2N (x/y), 



(16) 



where y = Ak + /(M N + E/A), and x = x Bj M N /M A , is 
increased by a factor of twenty towards the data as ex- 
trapolated in Ref . |24|] . But this calculation provides only 
a lower limit of the nucleon contribution because of the 
neglect of effects of the two-particle- two hole states |22| . 

Turn now to the experimental information about the 
nuclear pionic content. The Drell-Yan experiment on nu- 
clear targets pq ] showed no enhancement of nuclear pions 
within an error of about 5%-10% for their heaviest tar- 
get. No substantial pionic enhancement is found in (p,n) 
reactions [ p6[ . Understanding this result is an important 
challenge to the understanding of nuclear dynamics p7[ . 
Here we have a good description of nuclear dynamics, 
and our 5% enhancement is consistent [ps| , within er- 
rors, with the Drell-Yan data. 

This paper contains a new relativistic light-front the- 
ory of nuclear matter which leads to a good description 
of the binding energy, density, and incompressibility of 
nuclear matter. The use of a meson-nucleon Lagrangian 
enables us to also compute the mesonic content of the 
wave function using a consistent approach represented 
by Eqs. (0) , i) , (|) and @. The results are not in con- 
flict with extrapolations of deep inelastic scattering and 
Drell-Yan data to nuclear matter. 
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